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Abstract 

We study the coupled system of Maxwell and Dirac equations from a semiclassical 
point of view. A rigorous nonlinear WKB-analysis, locally in time, for solutions 
of (critical) order 0{y^) is performed, where the small semiclassical parameter e 
denotes the microscopic /macroscopic scale ratio. 

Key words: Dirac equation. Maxwell equations, nonlinear geometrical optics, 
WKB-asymptotics 

AMS (2000) subject classification: 81Q20, 35B25, 35B40, 35L60 



1. Introduction and Scaling 

The Maxwell-Dirac system (MD) is fundamental in the relativistic description of 
spin 1/2 particles. It represents the time-evolution of fast (relativistic) electrons 
and positrons within external and self- consistent generated electromagnetic fields: 

3 

ihds^ = a^{-ihcdk - e{Ak + *))^' + e{V + V''')^ + rnc^P^, 

(1.1) \ 

(\dl - a) y = Ip, f la? - a) A = ^j, 

J eo J ceo 

where the particle- and current- densities are defined by: 

(1.2) p--=em\ ife :=€(*, a'^vl-), fc=l,2,3. 

Here, ^I* = ^'(s, y) e is the 4- vector of the spinor field, normalized s.t. 

(1.3) / \^{s,y)\^dy=l, 

with s, y = (2/1,2/27^3), denoting the time - resp. spatial coordinates in Minkowski 
space. Further, F^'^^*) = V^"''^'' (s, y) e C is the self-consistent resp. external electric 
potential and A^*^^*^ = A^^^^\t,x) G C, the corresponding magnetic potential, with 
A = (Ai, A2, A3). In the following, the usual scalar-product for vectors X, F e 

^o-mail: christof.sparbor@univie.ac.at 
^e-mail: peter. markowich@univie.ac. at 



2 



will be denoted by {X,Y) and we shall also write \X\'^ :— {X,X). The so-called 
Dirac matrices /3, a'^, k = 1,2, 3, are explicitly given by: 



(1.4) /5 n _Tr I 



I2 \ k a'^ 
-I2/ ' " ■ I o-*^ 



with I2, the 2x2 identity matrix and a the 2x2 Pauli matrices, i.e. 

lA 2 -A 3 A 



Hence, a'^, are hermitian and moreover one easily checks that the following iden- 
tities hold for k,l — 1, 2, 3: 

f a'^a'- + a'- a'' = 24;, 
1 a''/? -t- (3a'' = 0. 



(1.6) 



Finally, the appearing physical constants are the normalized Planck constant h, the 
speed of light c, the permittivity of the vacuum eg, the particle mass m and the 
charge e. 

Additionally to Hl.l|l . we impose the Lorentz gauge condition 

(1.7) -^dsV{s,y)+divAis,y)^0, 

for the initial potentials V{0,y) and A{0,y). The gauge is henceforth conserved 
during the time-evolution. It ensures that the corresponding electromagnetic fields 
E, B are uniquely determined by 

(1.8) Eis,y):^--dsA{s,y)-VVis,y), B(s, y) curl A(s, y). 

c 

The MD equations are the underlying field equations of relativistic quantum electro- 
dynamics, cf. |Sc| . where one considers the system within the formalism of second 
quantization. Nevertheless, in order to obtain a deeper understanding for the inter- 
action of matter and radiation, there is a growing interest in the MD system also 
for classical fields, since one can expect at least qualitative results, cf. |EsSe| . 

From the mathematical point of view, the strongly nonlinear MD system poses a 
hard problem in the study of PDE's arising from quantum physics. Well posedness 
and existence of solutions on all of but only locally in time, has been proved 
almost forty years ago in |Gro| . On the other hand, only partial results (i.e. for 
small initial data) have been obtained in the quest of global- in-time solutions, |Ch| . 
[H] . |FSTj . let alone the study of other qualitative features of this system. 

In this paper, we shall analyze the MD system in a semiclassical regime. To do so, 
we first rewrite the equations such that there remains only one (positive) dimen- 
sionless parameter 

(1.9) ^ = ^, 

e^ 

which is obtained by replacing y y/y, s s/s and '^{s,y) y~^/'^'^{s/s,y/y), 
in order to maintain the normalization condition 11.3|l . with 

(1.10) cs — y and y — 



Here, we also replaced both, the external and the self-consistent potentials, by 
A(^^t){s,y) \A^'''=*'>{s/s,y/y) and F^-^^*) (s, y) ^ Ky('=^*)(s/s, y/y), with k = cX 
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and A = mc/e. We assume for the following that yj^e^;* yext g^,^ 0(1) in these 
units. In summary, we obtain the MD system in dimcnsionlcss form: 



(1.11) 



k=l 

UV = p, UA^ J, 



where from now on □ 9^ — A. In p.ll|l . s, y represent the microscopic time and 
length scales. Note that in general, S can not be considered as a small parameter, 
for example S « 13, in the case of electrons. Hence, semiclassical asymptotics in S, 
i.e. on Hl.ll|l directly, only make sense for highly charged and consequently heavy 
particles. 

Therefore, we need to rescale the system (|l.ll|l such that the time-evolution can 
be considered semiclassical, independent of the precise physical properties of the 
particles. We can suppose that the given external electromagnetic potentials are 
slowly varying w.r.t to the microscopic scales, i.e. V*^^* = V'^^'^{se/S, ye/ 5) and like- 
wise A*^^* = A'^'^*'{se/5, ye/ 5), where from now on e denotes the small semiclassical 
parameter. Here, the 5 is included in the scaling in order to eliminate it from the 
resulting equation. Hence, observing the evolution on macroscopic scales we are 
lead to: 

6 S 

(1.12) y^-x, s = -t. 

e e 

Moreover, we want that the coefficients of all nonlinearities to be 0(1), i.e. they 
should not carry a positive power of e. It turns out that there exists solutions , 
which obey this requirement. If we set 

(1.13) -■^{s,y)^,p'{t,x). 

e 

then, the normalization condition for (|1.3|l gives 

(1.14) / \r{t,x)\'dx^'- [ \^{s,y)\'dx^'-. 
This implies that we need to look for solutions -0^ s.t. 

(1.15) r-o{v^), 

assuming, as mentioned above, that d ^ 0{1). Wc therefore end up with the 
following semiclassical scaled MD system: 

3 

(1.16) iedtip' = a^i-iedk - (A| + A^^*))^'^ + /??/'^ + + V'^^W , 

fe=i 

(1.17) nv' = m\ 

(1.18) □A| = (V'^aV"), fc=l,2,3, 
subject to Cauchy initial data: 

V'1t.o = ^?(^)^o(Vi), 

(1.19) <( ^1t=o = ^/(^)' dtV^^^^Vfix), 

^^ = Af(x), dtA^l^^^AKx). 

For this nonlinear system, we want to find an asymptotic description of ^p^ ^ 0{^) 
as e ^ 0, i.e. a semiclassical description. Note that, equivalently, one could consider 
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asymptotic solutions of the form 

(1.20) ^O(l), 

which do not vanish in the hmit e ^ and which again satisfy the semiclassical 
scaled DM system, modified by the fact that the right hand sides of ()1.17|) . (|1.18|) 
are multiplied by an additional factor e. This illustrates the fact that we are dealing 
with a small coupling limit. We further stress that in our scaling the mass is 0(1) 
and fixed as e — > 0, which is different from the otherwise similar scaling used in 
|KS| . where a classical mechanics analogue of the DM system has been studied. 
The (rigorous) analysis of semiclassical asymptotics has a long tradition in quantum 
mechanics, the most common technique being the so called WKB-method. Quite 
recently, a semiclassical approach to the linear Dirac equation was taken in (BK) 
and also, using Wigner measures, in j(4MMP|, [Sp|. For a broader introduction on 
linear techniques and results, we refer to |SMMj and the references 

given therein. Nonlinear extensions of the WKB-method can be found for example 
in |Ge| . |Gr| . where scalar-valued semilinear Schrodinger equations are analyzed. 
We remark that the case of the nonlinear Dirac 4-system introduces significant new 
difficulties in the WKB-analysis, some of them are already present in the linear 
setting. 

Mathematically, our approach is inspired by techniques developed in |DoRa| , which 
sometimes go under the name weakly nonlinear geometrical optics. Due to the 
appearance of the small parameter e in front of each derivative in H1.16|l we are in 
the regime of so-called dispersive weakly nonlinear geometrical optics, which differs 
in several aspects from the non-dispersive one. We remark that the latter case is 
much better studied in the so far existing literature and we refer to I JMRSJ , for a 
recent review. 

To be more precise, we shall seek a local- in-time solution of (|1.16|l . which asymp- 
totically takes the following form: 

'4>^{t,x) = \/e u'^ {t, X, (j){t, x) / e), 
^^■^^^ ' u'{t,x,9)^ f^e^^^ujit^x^e). 

3=0 

Here, all the Uj{t, x, 9) g being 27r-periodic w.r.t. e M. Due to the factor e^/^, 
we call them small semiclassical approximate solutions, or small W KB- solutions. 

As expected, the phase function satisfies a (free) relativistic Hamilton- Jacobi 
equation and of course, convergence of the expansion H1.21|l . can only hold on a 
time interval, which corresponds to the existence of smooth solutions In weakly 
nonlinear geometrical optics, the homogeneity of the nonlinearity determines the 
required order of smallness of the asymptotic solution and, as we shall see, the 
factor ^/e precisely fits with the cubic nonlincaritics in Hl.l(j|l . We will show that 
for this particular scale we obtain, on the one hand, independent propagation of 
the electronic resp. positronic phase function ±0 and, the other hand, nonlinear 
interaction of the corresponding principal amplitudes uo.±- 

In other words, we study solutions on the threshold of adiabatic decoupling, a phe- 
nomena which is already well known in the linear case, cf. jPST| . In particular, the 
importance of the 0(^/e)-scale for the (linear) Dirac equation is stressed in |FKj . 
where one can also find a detailed description of the energy-transfer between elec- 
trons and positrons in terms of two-scale Wigner measures. These results, together 
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with ours suggest that if one wants to obtain semiclassical 0(l)-approximation in 
the strongly coupled regime, one needs to take into account simultaneously scales 
of order e and 0{^). These asymptotic solutions are then appropriate for heavily 
charged particles. We finally remark, that in the very recent paper |Je| . coupled 
Gauge-fields are studied from a similar point of view as in our work. 

This paper is organized as follows: We collect some preliminaries in section 2, then 
we shall determine the critical exponent and the corresponding eikonal equation of 
the approximate WKB-type solution in section 3. The corresponding e-oscillations 
introduced by the nonlinearity are determined in section 4 and the nonlinear trans- 
port of the approximation along the rays of geometrical optics is obtained in section 
5. Finally, in section 6 we shall prove that there exists a (local-in-time) solution of 
the MD system which stays close to the approximation and we also collect some 
further qualitative results. 

2. Preliminaries 

In the following, we will assume that no external electromagnetic fields are present: 



(2.1) y'=^*(t,a;) EE 0, A'=^*(i,j;) = 0. 
Moreover, we assume that at t = we have: 

(2.2) y/(a;) = V7(x) = 0, A]{x) ^ A]{x) = Q. 

Neither of these assumptions changes the following analysis significantly. They are 
only imposed for the sake of simplicity. Note that the DM system is time-reversible, 
but w.r.o.g. we shall consider positive times only in the sequel. 

Using the fundamental solution of the wave equation in dimension d = 3 and for 
times i > 0, we find the following expression for V , called the retarded potential: 

■^T^ J\x-y\<t \x-y\ 

= : gr{t,x)*\-^'{t,x)\'', 
where * denotes the convolution w.r.t. it, x) and 

(2.4) g^^t,x):=^6{t^\x\). 
Likewise, A can be written as: 

(2.5) A'[i^']{t,x) = Grit.x) * {r{t,x),o^r{t.x)). 

Using this representations, we shall rewrite (|1.16() - (|1.18|l in the form of a semilinear 
Dirac equation: 

j iedt^p" - V\{t, X, eD^)i}^ = 0, x e M^ i > 0, 

where I?^ is a matrix-valued differential operator [Dx := — iV). The corresponding 
e-dependent symbol is given by 

(2.7) V'j,{t,x,0^^-{i~A'W](t,x))+(3 + V'W](t,x), 
where x, ^, £ M"^, t e R. Here and in the following we use the notation 

3 

(2.8) a-i:=Y,a\u. 

k=l 



(2.6) 
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Note that in the nonhnear equation (|2.6(l . the potentials A'^[ip'^], V^[ip'^] depend 
non-locally on -0'^, as indicated by the bracket-notation. 

Multiplying H2.6|l with ip'^ and taking imaginary parts, we obtain the usual conser- 
vation law for a;)||2, hence the conservation of charge: 

(2.9) / ('0^(t,x),i/'^(t,x)) da; = |lV/(t,a;)|l2 = const. 

The free Dirac operator will be denoted by 

(2.10) V{eD^)iP^ -.^ -ieia-V)^" + P^P", 
with symbol 

(2.11) V{0=a-^ + (3. 

This 4x4 matrix has two different eigenvalues h±{(^) of multiplicity 2 each: 

(2.12) h±{0.^±m, ^eK', 

where 



(2.13) A(e) viei' + i, eeiR'- 

As expected, the eigenvalues h± (^) are nothing but the free Hamiltonian for a rela- 
tivistic particle. The positive resp. negative sign in H2.13|l corresponds to electrons 
resp. positrons. By straightforward calculations we obtain: 

Lemma 2.1. The spectral projectors n±(^) : C* C*, associated to h±{^) are 
given by 

(2.14) n±(0 i (^id4 ±^ v{0^ , n±n^ ^ o. 

The matrix- valued symbol I'(^) can therefore be decomposed into its positive and 
negative energy part in the following way: 

(2.15) vio = h+iou+io + h^m^io- 

Note that 

(2.16) h±{~0 = h±{0, whereas n±(-0 = U^iO- 
For later purpose, we also define: 

Definition 2.2. The partial inverse 

A±(0 : C ^ associated to U±{^) IS given 

by 

(2.17) A±(on±(0 = o, A±(oi?(C)x = (id4-n±(0)x, yxec\ 

Finally, we recall the definition of asymptotic equivalence: 

Definition 2.3. Let O C R", n > 1, be an open set, a^iy) e C°°(R";C") and 
g C'°°(]0, eQ[xIR"; C"). Then we say that a*^ is asymptotically equivalent to the 
formal sum X]j=o write 

oo 

(2.18) a^(y)^^£-'a,(y), 

3=0 

if for every m > 0, every multiindex a and every compact subset /C C O, there 
exists a C^.tr > 0, such that 



(2.19) sup 
/c 
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3. Generalized WKB-Ansatz and the Eikonal Equation 

At first we will show that the desired 0(Y^)-asyniptotics for the spinor field fits into 
the framework of weakly nonlinear (dispersive) geometrical optics, as introduced in 
|DoRaj ■ for nonlinear hyperbolic systems. 

We plug the following generalized WKB - Ansatz into equation l|2.()(l : 

' %l)^{t, x) = £^ u^t, X, (j){t, x)/e), 

^^'^^ ' w^(t,x,6») ~ ^e^fuj(t,a;,0), 

1=0 

where the functions Uj{t,x,9) E are assumed to be sufficiently smooth and 
27r-periodic w.r.t. e M. This gives 

(3.2) 

= ieP+^dtu'' - ePV%{t,x,eD^)u'' 
= ieP+^{dtu' + (a • W)u') - e^fiu' + ie^ {dt<t> + (a • y4>))dgu^ + e^PN^[u% 
with a nonlinear ity, J\f^ : — > C^, defined by 

(3.3) Af'iu"] {{a ■ A'[u']) - V^u"]) . 
The strategy is now to expand the right hand side of (|3.2|l as 

oo 

(3.4) eP^e^PRj{t,x) 

3=0 

and choose the coefficients Uj of (|3.1|l in such a way, that Rj{t, x) = 0, Vj G N. 

It is important to note that the first term on the right hand side of (|3.2() is of 
order eP+^, whereas the second and the third are ^ 0{eP). Since ^^[m^] 
are of order e^^*, by equation (|2.3() . (|2.5() . the function J\f'^[u'^] is of order e^^. This 
nonlinear term is supposed to be small, more precisely, it should not enter into the 
equation for Rq, describing terms of order 0(eP), but rather into expressions of 
0(e^+^). Thus we are led to the following normalization condition: 

(3.5) 3p = p+l, 

implying p = 1/2. With this normalization we have u"^ ~ 0(e^/^) (just as required 
by the scaling presented in the introduction), whereas the nonlinear term satisfies: 

AA^K] -0(£3/2). 

Remark 3.1. The choice p = 1/2 gives the critical exponent in the sense that 
for amplitudes 0(e^/^) one can prove simultaneously existence of the approximate 
smooth solution for times t = 0{1), i.e. on a time-scale independent of e, and 
nontrivial nonlinear behavior in the principal term of the approximation, cf. |DoRa| , 

Setting Ro{t, x) — 0, yields 

(3.6) i{dt(f>+{a-\/(l)))dgUo- (3uo = 0. 
Since Uj e C°°{M.'^ x §-'^;C'*) it can be Fourier- expanded w.r.t. 9 

(3.7) uo{t,x,9)^Y.''o,m{t^x)e'"''>. 

By this procedure, we find the following equation for the coefficients uo^rn'- 

(3.8) L{m d(j)(t, x))uo,m '■— {m dt4> + 'D{m'V(f))) mq.™ = 0, 
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where 'D{W{m(p)) is the 4x4 symbol matrix of the free Dirac operator evaluated 
at ^ = V(r7i(/)(i, a;)). In order to have a nontrivial solution uo,m 7^ we impose the 
condition, that there exists an open set fl C having a nontrivial intersection 

with {t = 0}, s.t. 

(3.9) det L(m d(j){t, x)) ^ 0, V(t, x) e n <Z R^+^. 
Using equations H2.12|l . (|2.15|l . this is equivalent to 

(3.10) imdtct)f -\mV(j)\'^ = 1, y{t, x) e n <Z R^+^ . 

Thus, for TO = ±1, the phase function (j) satisfies (in ft) the eikonal equation for the 
Klein-Gordon operator, i.e. 

(3.11) {dt(t)f - I V0p = 1, y{t, x)enc 

Indeed, it is easy to see that the choices to = ±1 are the only possibilities, since 
equation H3.11() gives 

(3.12) {mdt^f - |mV(/.p - 1 = (m^ - - 

= - 1, 

which is different from zero for all to 7^ ±1. Hence, in the Fourier-series (|3.7|) . there 
appear only two nontrivial harmonics, which are associated to the eikonal equation 
(|3.11|l : namely exp(i(/)(t,a;)/e), for to = 1 and exp(— i(/)(t, a;)/e), for to = —1. 
For TO = ±1 the equation (|3.1H) is fulfilled by two possible (p's, obtained from 



(3.13) 5t0±(t,x) = h±{W^±{t,x)) = ±V|V0±(t,x)|' + 1 = 0. 

This is the Hamilton- Jacobi equation for free relativistic particles. The following 
lemma guarantees existence and uniqueness of smooth solutions, where from now 
on, we shall denote by D"^ f{x), the Hessian of a given function / : R'^ M. 

Lemma 3.2. Given (t>i e C°°(M3;M), s.t. \\D'^4>i{x)\\ < C, there exist T± > and 
uniquely determined functions 4>± G C°°(r2-|-; M), where fi± := [0,T±) x R'^, s.t. 



(3.14) 



dt(b±{t,x) = h±{V(j)±{t,x)), y{t,x) e n±. 



Proof. We only proof the assertion for 0+, since the other case is completely anal- 
ogous. The initial value problem is non- characteristic everywhere, since 



(3.15) dt(l)+{0,x) = y^|V<?!>+(0,a;)|^ + l ^ 0, \fx e K^. 

Thus, dt(p+{Q, x) can be obtained from the initial data (t)+iO,x) = (j)i{x) G C^{R^) 
at each point a; e M^. Standard PDE theory then guarantees the existence of a 
unique smooth solution 0+ e C°°(il+;R), as long as 

(3.16) l-t\\D^h+{Vcj)i)\\ \\D^M ^0. 

Since D^H{^) is uniformly bounded, this condition holds by assumption and the 
assertion is proved. □ 

Remark 3.3. The assumption in lemma rOl can be relaxed to </>/ E C°°(M'^;M). 
In this case, however, one can not guarantee the existence a smooth solution 
in a space-time slab, but only in some open set O C R^"'"'^. In the following, this 
would lead to some technical difficulties, which we want to avoid, though, the whole 
procedure can be generalized to that case. 
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By equation H3.13|l . we have 

(3.17) -(/.±(t,x) = 

assuming that it holds initially a,t {t — 0} C fl. In the following, we therefore 
consider only the solution to (|3.13() with positive sign in front of the square root 
and write for it (j>{t, x) = (l)+{t, x). Also, we henceforth denote by := [0, T+) x M.^ 
the slab, in which existence of a smooth function (p is guaranteed. This we can do 
w.r.o.g. as will become clear in a moment: 

Since for <p = (j)^ it holds that dt4> — h^{Wcj)) = 0, equation (|3.8|l implies the 
following polarization conditions, locally for all (t, x) £ D,: 

(3.18) (n_(V0)uo,+i)(t,x) =0 ^ {Il+{\/(l))uo.+i){t,x)^uo,+i{t,x). 
Likewise, we get 

(3.19) (n+(V(/.)uo,-i)(t,a;) =0 ^ iU-{W(l))uo,-i){t,x) ^ uo,-i{t,x). 

One easily checks, using (|2.1t)|) and H3.17|l . that the conditions obtained with the 
choice = are equivalent to H3.18|l . H3.19|l . Thus, equation H3.8() indeed carries 
two degrees of freedom for the phase, given by ±(f> (or equivalently 0-). The 
amplitudes are then rigidly linked, by (|3.18|l . H3.19|l . 

In summary, we find that the principal term uo(t, x, 6) in our asymptotic description 
is given by 

(3.20) ua(t,x,(j){t,x)/e) := uo,+i(i, a;)e^'^(*'^)/^ + wo,-i(i, a;)e-'■'^(*•^)/^ 

where the amplitudes are polarized according to H3.18|l . (|3.19|) . From now on we 
shall use the simplified notation mo,±i = uo.±. 

4. Oscillations of the Nonlinearity 

Let us determine the response of the wave equations H1.17|l . H1.18|l to r.h.s. source 
terms induced by functions of the form H3.2()|l : 
To this end, we calculate: 

(4.1) |uo(i,x,(/.(i,a;)/£)|2 = \uoAi^x)? + |"o.-(^, a:)!'- 

The terms, which mix the electronic and positronic components cancel, since n± 
is hermitian and Ilillip = 0. Hence, we get from (12 .311 (at least formally), that the 
scalar potential V generated by the principal term mq, is simply given by 

(4.2) V[ua] = 5,(i,x)*(|uo,+ (t,x)|2 + |uo,-(t,x)p). 

In order to calculate the magnetic potential corresponding to uo, we first note that, 
by definition, we have the following identity 

(4.3) n±(e)(a-e + /3) = ^±(e)n±(c). 

Differentiating w.r.t. ^fc and multiplying (from the right) with n±(^) gives 

(4.4) n±(Oa'=n±(0 = n±(0(%/j±(0)n±(0 

= ± n±(0, 

since nj_(^) = n-i-(^). The expression 

(4.5) u:±{0:^V^h±{0^± ^ 
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is called the electronic rcsp. positronic group velocity, uj± e C°°(M'^; R'^). Using 
this definition, we obtain for k — 1,2, 3: 

(4.6) {uo{t, X, (t){t, x)/e), Q^uo{t, x, (j){t, x)/e)) 

= Lo+^k{y(^{t,x))\uo^+{t,x)\^ +Lo_^k{y(b{t,x))\uo,-{t,x)\^ 
+ {u„^+{t,x),a''uo,-{t,x))e-'^^^'^-^/' 
+ {uo^-{t, x), a''uo,+ it, a;))e''"*(*'")/^ 

The oscillating terms are usually called the Zitterbewegung of the Dirac-current, 
of. [Hfil, p. 195. The fact that the current-density corresponding to uq carries 
e-oscillations is in sharp contrast to the WKB-approach for Schrodinger-type prob- 
lems, see e.g. |Gej . |Gr| . 

The Zitterbewgung may cause severe problems since a-priori one cannot exclude the 
possibility of resonant interactions between the principal term uq and the magnetic 
potential A^[mo] obtained from H1.18|l with r.h.s. given by (|4.6|) . If this happens to 
be the case, our one-phase ansatz H3.1(l breaks down and instead one would need 
to establish a so-called resonant asymptotic expansion in the spirit of I.TMRlj . (We 
remark that so far, only the case of resonances in one spatial dimension can be 
treated rigorously, cf. |JMR2| .') 

We will show that these problems do not appear in our situation. To this end, we 
need to describe precisely what kind of e-oscillations are present in [uq] . 

First we note that, by the superposition principle, every term appearing on the 
r.h.s. of 1)4. 6|l generates its own potential field. The nonoscillating terms of (|4.6|) 
lead to a standard hyperbolic problem, hence H2.5|l gives 

(4.7) Aa[uo]{t,x) := ^^.(t,^) * (cc;+(V(/))|mo,+ |' + c^-,fe(V0)|uo^_|') {t,x). 
In order to treat the Zitterbewegung, let us define Z :— [Zi, Z2, Zj,) by 

(4.8) Zk{t,x,(j,{t,x)/e) := (wo,+(i, a^), a'"o,-(t, a;))e-'20(M)/e 

+ (uo,-(t,a;),a'=uo,+ (i,2:))e^20(M)/e^ fc= 1,2,3. 
Using this definition, we can now prove the following lemma: 

Lemma 4.1. Let f2 C be the slab in which existence of a smooth phase 

(j) e C°°{n;R), satisfying WW , is guaranteed. Then, given Z € C°°(r2 x S^C^), 
as in \4--<^ , there exists a uniquely determined smooth A'^ e C°°{fl x §^;C^), with 

00 

(4.9) A^t,x,e) r^J2^^Mt,x,s), 

1=1 



s.t. A^{t,x,(j){t,x)/e) satisfies: 

( aA^t, x) - Z{t, x) - 0, in C°°(rJ; C^), 
1 A^\^d,A^\^ 0. 



(4.10) 



it=o 



More precisely, Ai{t,x,<f>{t,x) / e) can he written in the following form: 

(4.11) A,(t,a;,(/.(i,a;)/£) = A+(t,x)e'20(t,.)/e^^-(^^^)g-^20(t,.)/e^ 

with principal amplitudes A^ £ C°°(r2;C^), given by 

(4.12) A^.ki^,x) ^ ± {uo{t,x),o^uo{t,x)). 
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Proof. The proof can be done separately for each spatial component of A^{t,x,9) 
and for both types of oscillations, corresponding to ±.2(j). Hence, we are lead to the 
following type of problem: 

I aa'it, x) = h{t, (t, x) G 

l«1t=o=^*«1t=o= 0, 
for some given h E C°°{^; C). Let us define a new variable f^{t, x) G by 

f{t,x) := (9la^,92a'',93a^,(9ta^a^)^(^,a::). 
Further, denoting by 

b{t,x) := (0,0,0,6(t,a;),0)^, 
we can rewrite (I4.13f) in the form of a symmetric hyperbolic system 

(4.14) {dt + iX-V) + K) fit, x) - b{t, :r)e±^20(t,.)/s^ 

with A'^, k, denoting realvalued (symmetric) 5x5 matrices. In our case, these 
matrices are simply given by (see e.g. p , p. 21, for more details): 

(4.15) := - ((5mfc(5„4 + 5„i4<5nfc)™^„ , 1^ ■= - {5^m5An)rn.n^ m, 71 = 1, . . . , 5, 

where 5ab denotes the Kronecker symbol and k = 1,2,3. 

It is now possible to use the existing results on linear geometrical optics, provided 
the phase ±20 is not characteristic for the system (|4.14|) . i.e. 

(4.16) det(±2at0(t,a:) ±2A- V0(i,a;) + k) 7^ 0, V(t,a;)efJ. 
Computing this determinant, we obtain the condition 

(4.17) ±32(9t</.)3((5t0)2- iV^n ^0, V(i,a;)ef]. 

Since, by assumption, (j) solves the Klein-Gordon eikonal equation H3.11|l the second 
factor on the l.h.s. of (|4.17|l is equal to one and thus, different from zero in all of 
fi. On the other hand we get from (|TTT|) : {dt(j)f = (| V^ip + l)^/^ ^ 0, V(t, x) G n. 
Hence, condition (|4.16|) is fulfilled and the assertion follows from theorem 4.4 in 
|H,a| . In particular we get: 

(4.18) Af^kii^x) = _^Q^^y\ |y^|2 {uo{t.x),a''uo{t,x)) 

= ± (uo(t, a;),a''Mo(t, a;)), 
which concludes the proof. □ 

Lemma [4.11 shows that the Zitterbewegung in (|4.6() generates a magnetic potential 
which is small, i.e. at least of order 0{e). Moreover the warepszZon-oscillations, 
appearing in A'^ , are exactly the same as in the (|4.6|) and hence we can consistently 
proceed with our one-phase expansion for '0'^. 

Remark 4.2. Although the MD system is hyperbolic, lemma HITl can be considered 
as an analogue of so-called elliptic high frequency asymptotics, the main feature of 
which is the fact that asymptotic solutions can be obtained by local (in t, x) alge- 
braic relations. In other words, the Maxwell system can be considered transparent 
w.r.t to the oscillations generated by the Dirac equation. 

The result of lemma 14.11 implies that the nonlinearity [uq] , defined in H3.3|l , 
admits an asymptotic expansion of the form 

oo 

(4.19) U'[uo]{t,x,6) ^ Mo[uo]{t,x,6) + Y.e'Mi{t,x,6), 

1=1 
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where, using the expressions (|4.2(l . H4.7|l . we have: 

(4.20) no{t,x,e)= {{a- Ao[uQ]{t,x))-V[uQ]{t,x))uo{t,x,e) 
and for alH > 1: 

(4.21) Ni{t,x,e):= [a- Mt,x,e))uo{t,x,e), 
with Ai given by (|4.11|) . 

Note that the expression H4.19|l represents two kinds of e-oscihations: Those de- 
scribed by phase-functions ±0 are present in aU terms Afi with I > 0, whereas 
e-oscihations with phases ±3(/) appear in Afi with I > 0. Also, note that uo enters 
in a nonlocal way only in the lowest order term H4.2n(l . 



5. Nonlinear Transport along Rays 

We need to find an evolution equation, which determines uq from the initial data. 
To this end, let us define an operator P, which projects on the set of harmonics 
corresponding to solutions of the eikonal equation 13.11|l : 

Definition 5.1. Given some v G C°°(M^ x §1;^), which can be represented by 

(5.1) i;(i,a;,e) - ^i;„(i,a;)e™^ 

we define the action of P on u, by 

(5.2) {Vv){t,x,e) := (n+(V0)w+i)(i,x)e^« + (n_ (V(/))i;_i)(t, x)e-'^ 

In words: P picks modes corresponding to m = ±1 and multiplies them with the 
matrices n±(V0). Note that, at least in H., it holds true that 

(5.3) i¥uo)it,x,9)^uo{t,x,e), 
in view of (EiHll and (pTT^ . fTTi^ . 

From (|3.2() . we have that the evolution of uq is determined by terms of order 
0(eP+^) — 0{e^^^). Setting the corresponding coefficient in H3.4|l equal to zero, 
i.e. R2{t,x) — 0, yields 

(5.4) i{dt(t) + (a ■ V<l)))deU2 - (iu2 + iidtuo + (a ■ V)uq) + 7Vo[uo] = 0, 
with A/'o[uo] as in H4.20|l . Equation H5.4|l imphes 

(5.5) iidt + (a ■ V))uo+Ma[uo] £ ran (z(9t0 + (a ■ Vcj}))de - /?) . 

Applying P to H5.4|l . eliminates the term including U2, since P projects on the kernel 
of i{dt<j) + {oi ■ V(f)))d0 — /3 and we obtain 

(5.6) iWtUQ + if {a ■ V)uo + P7Vo[mo] = 0. 
Using the fact that Puq — uq^ by l|5.3(l . this gives 

(5.7) P9t(Puo) + P(a • V)(Pmo) = iPA/'o[Puo]- 

This equation is similar to the one appearing in jPoRaj . however, in contrast to 
the quoted work, our nonlinearity constitutes only the first term of an asymptotic 
expansion of the full [uq] ■ 

We proceed by stating a useful identity: 

(5.8) a'n±(0 = n^(0«' + ^±,fc(0l4, 
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obtained from straightforward calculations. After more lengthy but straightforward 
calculations, in which we apply the relations H4.3|l . (|4.4(l and H5.8|l . we can express 
the l.h.s. of H5.7|) in the form of a transport operator: 

(5.9) n±at(n±uo) + n±a- v(n±uo) = 

9tUo,± + (w±(V(/)) • V)mo,± + ^ div{uj±{W(j)))uo,±. 

On the other hand, computing the action of the projector P on the nonlinear term 
7Vo[uo], we get 

(5.10) PAAopuo] - e'*/' {{Ao[uo] ■ a;+(V0)) - V[uo]) uo.+ 

+ e-'^/^ {{A„[uo] ■ c^-(V0)) - V[uo]) uo,- . 

Here we have again used (|4.4() . Thus, we finally conclude, that the time-evolution 
of the principal amplitudes uo,± is governed by the following semilinear first-order 
system: 

^^■^^^ j(at + (c^-(V(/))-V))uo^_(i,.x)= T_[uo]{t,x)uoMt,x), 

where 

(5.12) r±[uo]{t,x) :=^^oK]•c^±(V0)-^F[uo]-^div(cJ±(V(^)). 

By construction, the polarization of uq ± is conserved during the evolution. The 
system (IS.llf) determines (Puo)(t, x, 9), from its initial data (Puo)(0, x, 9) and since 
(Puq) = uoj we have completely constructed uq. 

Multiplying (jOT)) by uo,+ resp. uo,~ and integrating by parts, we obtain the 
important property of charge-conservation: 



(5.13) / \uQ^-f-(t, x)\ +\uo^-{t,x)\ da; = ||uo(i, a;)||2 = const. 

Given uq, determined by H5.11|l . it remains to construct the higher order terms 
Uj(t,x,9), J > 1 of our approximate solution. This can be done by a similar 
construction as given in jPoRaj : 

We expand the cubic nonlinearity M^[u^] in powers of e: 

(5.14) A/'^[-ito + H ]^ W[uo] + eM^[uo,ui] -\ , 

where, using the definitions l|4.2|l . 1)4. 8|l . we easily compute: 

(5.15) M'^[uq, ui] = 2 {Qr * {uo,ui))uq + V[uo]ui + a ■ {Gr * Z) ui 

We need to apply lemma UTI to all terms appearing on the r.h.s of H5.14|l . which 
results in a similar expansion as given in H4.19|l . Hence, after rearranging terms in 
powers of £, we can write 

(5.16) 7V=[uo + eui-| ]- WoK] + eA^i + eA^oK, "i] + • • • ■ 

Consequently, for j > 1, the 0(e-'/^+^/^)-coefiicient is given by 

(5.17) Rj{t,x) = i{dt(t)+ (a ■ "^4>))deUj - Puj + i{dt + {a ■ V))uj_2 

+ Mo[uo,Uj-2] + S{ua, . . . ,-u„<j-2), 
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where, as usual, we impose: Un{t,x,6) — 0, for all n < 0. The source term 5, 
only depends on lower order coefficients uq, . . . , Un<j-2- It is obtained by applying 
lemma I^TI to higher order terms in the expansion (|5.14|l . leading to contributions 
Ml and Ml with: / + 1 = j/2. 

We can now decompose 

(5.18) u^{t,x,e) = {Puj){t,x,9) + {idi-F)uj{t,x,9). 
Note that in contrast to uq, where, in view of H5.3|l . it holds 

(5.19) {id4-P)uoit,x,e) :^0, 

we can not expect all higher order coefficients Uj to be polarized too. Hence, we 
need to determine separately Puj and (id4— P)itj. To this end, we introduce the 
following definition: 

Definition 5.2. Again, let v{t, x, 9) be given as in definition l5.1l then we define a 
partial inverse Q, associated to P, by 

(5.20) {Qv){t,x,9) := (A+(V0)«+i)(i, a;)e^^ + (A_ (V0)w_i)(t, a;)e-*^ 
where A± is the partial inverse to n±, defined by (j?T7|) . 

Assume now that we already know m„, for n < j, then (id4 —P)uj is determined by 
setting {QRj){t,x,9) = 0. This gives 

(5.21) 

(id4 -P)uj = - Q(^i{dt + (a • V))uj_2 + A^oK, "^-2] + S{uq, . . . , m„<j-2)) • 

On the other hand, setting {PRj^2){t,x,9) = 0, we obtain an evolution equation 
for Puj-. 

(5.22) iPdtiPuj) + iP{a ■ VPwj) = ~PMo[uo, u^] + r(uo, . . . , uj), 
where 

(5.23) r(uo, . . . , Uj) := -PS{uo, . . . , m„<j) - P + ^(a • V) - /?) (id4 -P)uj. 

Here, the first term on the r.h.s is already known by the inductive hypothesis and 
the second one is given by equation (|5.21|1 . Hence, by induction, one can construct 
all higher order coefficients Uj{t,x,9), j >1 in this way. 

Note, that the left hand side of H5.22|l is essentially a transport operator, which can 
be expressed as shown above. Thus (I5.22|l constitutes a linear first order system, 
which determines the so-called propagating part Puj from its initial data. 

Remark 5.3. The above construction can be generalized to the case, where, ad- 
ditionally given external potentials A*^^', are included, or, equivalently, non- 
zero Cauchy initial data for the Maxwell equations 11.17|l . 11.18|l are assumed. One 
checks that, instead of H3.13|l . the following Hamilton- Jacobi equation, correspond- 
ing to m = 1, holds: 

(5.24) 9t(/)± ± ^J \\/(f>± - A^^*{t, x)\^ + 1 + V^'Ht, x) = 0. 

Since no other harmonics with to 7^ 1 exist, one again ends up with two phases 
(t)±{t,x), corresponding to the electronic resp. positronic degrees of freedom. In 
this case however, —(j)+{t,x) ^ (j)-{t,x), in contrast to (|3.17|l . Also, one obtains 
an additional matrix-valued spin-transport term, appearing on the left hand side of 
(|EjT|| and which can be found in [EKl, EH, e.g.. 
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We are now in the position to formulate our first theorem (in which we do not aim 
to impose the weakest possible assumptions). In the following, C"^^ denotes the 
space of smooth function, compactly supported in a; G M"^. 

Theorem 5.4. Assume that the initial data ^pi{x) admits an asymptotic expansion 
of the form: 

ipj{x) = y/eu''{x,(j)i{x)/e), 



(5.25) 



u^(x,0)~^e^/2x,(.T,0), 



(5.28) 



where (Pi G C°°(R3;]R) satisfies \\D^(Pi\\ < C. Further, let Xj e C^){^^ x S^C^) 
be s.t. 

(5.26) {¥xj)i^,0)=xjix,e), VjeN. 

Then, there exists aO <T* <T, a corresponding domain fl* :~ [0,T*) H Q and a 
uniquely determined u'^ e C'^^{fl* x §^;C'*), with 

oo 

(5.27) u^{t,x,9) - ^/e^e^''^Uj{t,x,e), 

3=0 

s.t. u'^[t,x,<j){t,x)/e) satisfies: 

j isdtu^ -V%{t,x,eD)u^ - 0, V(t, x) e n* , 

More precisely we have: 

The principal term uq is given by ^MW, satisfies (Puo)(<, x,9) = uo{t,x,9) and 
solves j5.7| ) with initial data (Puo)(0, x, 0) = xo{x,9). 

For all j > 13, the infinite sequence of equations \5.21]) . H5.2il\} . uniquely determines 
Uj{t,x,9), with initial data {Puj){0,x,9) = Xj{x,9). 

Proof. The existence of a smooth phase G C°°(fi;M), on the slab C Ri+'^, is 
already guaranteed by lemma 

Next, consider the case j = 0: Since ijJ±(V0) G K^, defined by (|4.5|l . satisfies for all 
multiindices cr, v: 

(5.29) sup \d^d^uj±,kC^(l){t,x))\ <oo, /c = l,2,3, 

we find that the l.h.s. of (|5.11|) constitutes a linear symmetric hyperbolic system. 
From i^-conservation property H5.13|l the usual commutator estimates lead to H'^- 
regularity, i.e. uo G C^{M.'^; H'^) for all s > 0. Now, it is a standard results for 
the linear wave equations in d = 3 spatial dimensions, that source terms in iJ*(M'^) 
generate solutions in H^{R^), cf. ,H0) , chapter XXIII. This fact and Schauder's 
lemma imply that the maps 

(5.30) uoit,-)^r±[uoit,-)] 

are locally Lipschitz from H'^{{0,t) x x S^) to itself, for all s > 2, uniformly 
for < t < T. By a standard Picard iteration we therefore obtain a local-in-time 
existence and uniqueness result in H''{V,* x §^), for every s > 2 and a Sobolev 
imbedding gives uq G C^{il* x S^; C*). The proof of the asserted regularity for the 
i-derivatives follows by using the differential equation to express them in terms of 
x-derivatives and the finite speed of propagation for solution of ()5.11f) implies that 
Uq is compactly supported in Rj^ since xo is. 
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Finally, for j > 0, we have that the amplitudes Uj^± are determined by the linear 
symmetric hyperbolic system H5.22|l . (15 .2111 and the assertion is proved. □ 

Once again, we stress the fact that we analyze the MD system in a weakly coupled 
regime. Indeed, the above result implies: 

Corollary 5.5. Let u'^{t,x,(l){t,x)/e) be as in theorem \5.4\ then 
( V'[u']{t,x) ^ eV[uo]it,x) +0{e^), 
\ A' [u'] {t, x) - eAa [uo] [t, x)+Oie^), 
where V[uq], ^o[^^o] nonoscillating and explicitly given by ^4- 



(5.31) 



6. Stability and further results 

In theorem 15.41 we obtained a function u^, which solves the MD equation up to a 
residual i?^ ^ 0, compactly supported in [0, T*] x R^. We want to compare to a 
true solution ip^ and prove that u'^ — ^ on Q* ^ [0, T*] x R'^. 

Theorem 6.1. Under the assumptions of theorem \5.4\ there is an e* E (0, 1), s.t. 

for e < s* , there exists a unique smooth -0^ G C^^^(ri*; C*), satisfying 

j iedtip" -V%{t,x,eD)iJj'' = 0, y{t,x) e fl* , 



(6.1) 



which is asymptotically equivalent to u^, i.e. 

(6.2) ■iP%t,x) u%t,x,(l3{t,x)/e) m C(^)(f7*;C4). 

Proof. Defining v"^ := u'^ — , we obtain for the following IVP: 

i{edt^e{a-V))v^ - Pv^ ^N'[u' + v^]~N'[u']= - R\ in f7* 

0, 



(6.3) 



it=o 



The nonlinearity can be handled analogous to the proof of lemma 6.2, in |DoRaj . 
since for smooth sources the wave equation has smooth solutions, which moreover 
travel with finite speed. Having this in mind, the rest of the proof is a simple 
modification of the one of theorem 6.1 in jDoRaj . □ 

As far as particle creation is concerned, the local-in-time solution t/j^ ^ 0{-\/£) 
shows the following qualitative behavior: 



Corollary 6.2. Let V/ be as m theorem\5^\ If (n-ipj){x) = 0, then for < t < T* 
it holds: (n_i/'^)(i,a;) ~ 0{e^/'^), i.e. no positrons are generated, up to 0(e^/^) and 
the analogous statement for electrons is valid, too. 

Proof. The assertion holds true, since a careful examination of the asymptotic ex- 
pansion shows that both, uo and wi, satisfy: {¥uj){t,x,6) — Uj{t,x,9), in il*. □ 

For completeness we shall also consider the matrix-valued Wigner transform corre- 
sponding to ip^, i.e. 

1 f f . /. e 

(2^)3 

where ® denotes the tensor product of vectors. The 4 x 4-matrix is a phase- 

space description of the quantum state ip"^ . 



(6.4) m (t, x,0-= 7^ li;'(t,x+^-y)®r- (t, X - |y) e^^'^dy. 
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Corollary 6.3. Let '0'^ ^ 0{^) be the unique smooth local-time- solution of the 
MD system, as guaranteed by theorem j6'.-?|) and let w^lip'^] ~ 0{e) be its Wigner 
transform. Then, up to extraction of subsequences, we have 

(6.5) lim '^w^l'iij^] = n in S'{[0, T*) x x R^) weak-k, 

where the matrix-valued Wigner measure fj, is given by: ji = /i+ + with 

(6.6) fj.±{t, X, C) = uo.±{t, x) «) uo,±{t, x) 6{(, =F V0(i, x)). 

Proof. Since (/) has no stationary points within Q* , a nonstationary phase argument 
imphes that all Wigner matrix elements, which mix the electronic and positronic 
components are of order 0(e°°). The assertion then follows from the well known 
results on Wigner measures, cf. iGMMP . □ 

We finally remark on the case of the Dirac-Maxwell system where the Dirac particles 
have vanishing mass. Instead of (|3.13l) we obtain 

(6.7) a*^ = ±|V0|, 

which is equivalent to the eikonal equation of the wave equation. It follows that 
in this case lemma I^TTl can not hold, since the phases ±(j> are characteristic for the 
wave equation. More precisely, they are indeed everywhere characteristic, i.e. in all 
of fi, which again allows for an asymptotic description of the A^, similar to H4.9|l . 
I)4.11|l . cf. jLa| or |Raj . chapter 5. In this case, the e-oscillations are also given by 
exp(±2i(^/e), but the corresponding amplitudes Af" are of course different. The 
main difference, however, is the fact that in this case the summation index runs 
from I = to infinity, i.e. e-oscillations are present already in the lowest order 
term. This leads to a more complicated structure of the transport equations for 
the amplitudes uj, but apart from that all results remain valid. 



Acknowledgement : 

This work originated from a discussion with P. Gerard (Univ. do Paris Sud). The 
authors thank J. Ranch and H. Spohn for helpful discussions and the TICAM 
(Univ. Texas) for its support. This work has also been supported by the Austrian 
Science Foundation FWF through grant no. W8, project no. P14876-No4 and the 
Wittgenstein Award 2000 of P. M. Additional financial sponsorship has been given 
by the European Union research network HYKE. 

References 

[BK] J. Bolte, S. Keppeler, A semiclassical approach to the Dirac equation, Ann. Phys. 274 
(1999), 125-162. 

[Ch] J. M. Chadam, Global solutions of the Cauchy problem for the (classical) coupled 

Maxwell-Dirac equations in one space dimension, J. Fund. Anal. 13 (1973), 173-184. 

[DoRa] P. Donat, J. Ranch Dispersive nonlinear geometrical optics, J. Math. Phys. 38 (1997), 
1484-1523. 

[EsSc] M. Esteban, E. Sere, An Overview on linear and nonlinear Dirac equations, Discrete 

Contin. Dyn. Syst. 8 (2002), no.2, 381-397. 
[FK] C. Fermanian-Kammerer, Semi-classical analysis of a Dirac equation without adiabatic 

decoupling, preprint, to appear in Monatsh. f. Math. (2003). 
[FST] M. Flato, J. C. H. Simon, E. Tafflin, Asymptotic Completeness, Global Existence and 

the Infrared Problem for the Maxwell-Dirac Equation, Memoirs of the AMS no. 127. 
[G] V. Georgiev, Small amplitude solutions of the Maxwell-Dirac equations, Indiana Univ. 

Math. J. 40 (1991), no. 3, 845-883. 



18 



[Ge] P. Gerard, Remarques sur I'analyse semi-classique de I'equation Schrddinger non 

lineaire, Seminaire EDP de L'Ecole Polytechnique (1992), lecture no. XIII. 
[GMMP] P. Gerard, P. Markowich, N. Mauser, F. Poupaud, Homogenisation Limits and Wigner 

transforms, Comm. Pure Appl. Math. 50 (1997), 323-379. 
[Gr] E. Grenier, Semiclassical limit of the nonlinear Schrodinger equation in small time, 

Proc. AMS 126 (1998), no.2, 523-530. 
[Gro] L. Gross, The Cauchy problem for the coupled Maxwell and Dirac equations. Comm. 

Pure Appl. Math. 19 (1966), 1-15. 
[Ho] L. Hormander, The Analysis of Linear Partial Differential Operators III, Springer 

(1985). 

[Je] P. Y. Jeanne, Optique geometrique pour des systemes semi-lineaires avec invariance de 

jauge, Mem. Soc. Math. Fr. (N.S.) 90 (2002). 
[JMRl] J. Joly, G. Metivier, J. Rauch, Resonant one dimensional geometric optics, J. Funct. 

Anal. 114 (1993), 106-231. 
[JMR2] J. Joly, G. Metivier, J. Rauch, Recent results in nonlinear geometric optics. Internal. 

Series Num. Math. 130, Birkhauser (1999). 
[KS] S. Kunzo, H. Spohn, Adiabatic limit for the Maxwell- Lorentz equations, Ann. Henri 

Poincare 1 (2000), no. 4, 625-653. 
[La] P. Lax, Asymptotic solutions of oscillatory initial value problems, Duke Math. J. 24 

(1957), 627-646. 

[MaFe] V. P. Maslov, M. V. Feydoriuk, Semi- Classical Approximation in Quantum Mechanics, 

Reidel Dordrecht (1981). 
[PST] J. Panati, H. Sohn, S. Teufel, Space- adiabatic perturbation theory, preprint 

math-ph/0201055 available at: |http://xxx. lanl.gov 
[R] R. Racke, Lectures on Nonlinear Evolution Equations, Vieweg (1992). 

[Ra] J. Rauch, Lectures on Nonlinear Geometrical Optics, IAS/Park City Math. Series 5, 

AMS (1999). 

[Ro] D. Robert, Semi-classical approximation in quantum mechanics. A survey of old and 

recent mathematical results, Helv. Phys. Acta 71 (1998), 44-116. 
[Sc] F. Schwabl, Advanced Quantum Mechanics, Springer (1999). 

[SMM] C. Sparber, P. Markowich, N. Mauser, Wigner functions vs. WKB-methods in multi- 
valued geometrical optics, Asympt. Anal. 33 (2003), no.2, 153-187. 

[Sp] H. Spohn, Semiclassical limit of the Dirac equation and spin precession, Ann. Physics 

282/2 (2000), 420-431. 



